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The adiabatic connection of density functional theory (DFT) for electronic systems is generalized 
here to negative values of the coupling strength a (with attractive electrons). In the extreme limit 
a — > — oo a simple physical solution is presented and its implications for DFT (as well as its 
limitations) are discussed. For two-electron systems (a case in which the present solution can be 
calculated exactly), we find that an interpolation between the limit a — > — oo and the opposite 
limit of infinitely strong repulsion (a — > +oo) yields a rather accurate estimate of the second-order 
correlation energy E^ h2 [p] for several different densities p, without using virtual orbitals. The same 
procedure is also applied to the Be isoelectronic series, analyzing the effects of near-degeneracy. 



I. INTRODUCTION, DEFINITIONS AND 
OUTLINE 

Combining low computational cost with reasonable ac- 
curacy for many molecules and solids, density functional 
theory (DFT) [1, 2] has become a particularly success- 
ful approach for electronic-structure calculations, both 
in chemistry and in physics. In DFT, the exact electron 
density p = p(r) and ground-state energy En[v] of N in- 
teracting electrons in a given external potential v = v(r) 
could be in principle obtained by solving single-particle 
("Kohn-Sham") equations for non-interacting electrons. 
In the practical implementation of Kohn-Sham (KS) 
DFT, however, we have to rely on approximations for 
the density functional i? xc [p] of the exchange-correlation 
(xc) energy. Despite the large number of available ap- 
proximations for this functional and of their success- 
ful applications, there are still important cases in which 
KS DFT can fail, which is why the quest for better xc 
functionals continues to be a very active research field 
(for recent reviews see, e.g., Refs. 3-5). For example, 
present-days KS DFT encounters problems in the treat- 
ment of near-degeneracy effects (rearrangement of elec- 
trons within partially filled levels, important for describ- 
ing bond dissociation but also equilibrium geometries, 
particularly for systems with d and / unsaturated shells) , 
in the description of van der Waals long-range interac- 
tions (relevant, for example, for biomolecules and layered 
materials) , and of localization effects due to strong elec- 
tronic correlations, as those occurring in Mott insulators 
or in low-density nanodevices. 

An exact expression for the xc functional is the 
coupling-constant integral [6-8], 

E xc [p] = f daW a [p}. (1) 
Jo 

The integrand is the a-dependent density functional 



with the operator of the electronic Coulomb repulsion, 




and the continuum functional of the Hartree energy, 

The crucial quantity in Eq. (2) is the a-dependent wave 
function 

*<*[p] =V a {[p];r 1 ,...,r N ;a 1 ,...,a N ), (5) 

where the r„ and the cr„, respectively, are spatial and spin 
coordinates of the electrons. Out of all antisymmetric 
A^-electron wave functions 'F that are associated with the 
same given electron density p, *& a [p] denotes the one that 
yields the minimum expectation of T+aV ee [9], with the 
kinetic-energy operator T = —^mT SilLi , 

(y a [p}\f + aV ee \y a [p]) = min(*|f + aV cc \V). (6) 

If the density p is w-representable for all a > 0, there 
exists an a-dependent external potential v" xt ([/?],r) such 
that ^f a [p] is the ground state of the Hamiltonian 

N 

H a [p] =f + aV cc + J2 <Mri). (7) 

i=l 

By construction, this Hamiltonian has for all a > the 
same ground-state density p = p(r) as the real system 
with a = 1. 

At the non-interacting limit a = 0, the ground state 
of H a [p] is, in most cases, the single Slater determinant 
^o[p] with the N occupied Kohn-Sham orbitals. Conse- 
quently, 



W a \p] = (* a \p]\V ee \* a \p])-U\p], 



(2) 



E x [p] = W [p] = (^o[p}\V cc \^ [p}) - U[p] (8) 
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is the functional of the DFT exchange energy. 

For a > 0, 'fcj/o] is no longer a Slater determinant of 
single-particle orbitals, but a correlated TV-electron wave 
function. As a > grows, the electron-electron repul- 
sion in the state ^ a [p] increases. So does the average 
distance (|rj — rj\) between two electrons. Consequently, 



for N > 1, the expectation of V cc in Eq. (2) must be a 
motonically d< 
The quantity 



monotonically decreasing function of a, j^W a [p] < 0. 



1 d 



a=0 



(9) 



is the 2nd-order correlation energy in the Gorling-Levy 
perturbation expansion [10, 11]. It can be expressed in 
terms of the KS single-particle orbitals, but, in contrast 
to E^lp], it requires also all the unoccupied orbitals, 



\(^[p]\V ee -Vu[p]-V^[p} \MP})\' 



k=i 



Eh - En 



(10) 

Here, ^q[p] is the fc-th excited state and Eq the corre- 
sponding eigenvalue (while En is the ground-state energy) 
of the non-interacting Hamiltonian Hn [p] . The operators 
Vn[p] = E^iM[p];ri) and V x [p] = ^ 1 v x ([p\; r;), re- 
spectively, represent the Hartree potential, 



M[p\;r) 



SU[p] 



5p(r) 

and the exchange potential, 
*>x([p];r) = 



2 f ^ pOO 



5EM 
Sp(r) ■ 



(11) 



(12) 



Since the exchange functional E x [p] is not known explic- 
itly in terms of the density p, but only implicitly via 
the KS orbitals [Eq. (8)], the evaluation of the function 
w x ([p];r) for a given density p is a non-trivial problem 
(see, e.g., [12-17]). The resulting weak-interaction ex- 
pansion of W a [p] is then 

W a [p] = E x [p] + 2E? L2 [p]a+... (a-0). (13) 

The functional VF a [yo] and its exact properties have al- 
ways played a central role for the construction of approx- 
imate E xc [p] (see, e.g., [18-22]). Although the integra- 
tion over a in Eq. (1) runs between and 1, we can 
consider values of a larger than the physical interaction 
strength, a > 1. In particular, the strong-interaction (or 
low-density) limit of DFT is defined as the a — > oo limit 
of W a [p] . It has been shown that in this limit the leading 
terms in W a [p] are [23-27], 



W \p] W" \p\ 



a 



(a-»oo). (14) 



While generally W^[p} = [23, 27, 28], the coefficients 
Woo[/9] and W^Jp] have been evaluated systematically 



for spherically symmetric A-electron densities [26, 27]. 
This a — > oo expansion of VF Q [/9] is useful in several 
ways. For example (even if treated in an approximate 
way [23]), it has been used to build an exchange corre- 
lation functional (interaction-strength interpolation, or 
ISI) by interpolating W a [p] between a — > [Eq. (13)] 
and a — > oo, yielding atomization energies of simple 
molecules with errors within ~ 3.4 kcal/mol [29, 30]. 
Moreover, by properly rescaling a given approximate 
functional E xc [p], it is possible to test its performance 
in the strong-interaction limit, thus adding a new con- 
straint for building approximations [28, 31, 32]. More 
recently, the strong-interaction limit has been used to di- 
rectly address strongly correlated systems [33, 34] in a 
DFT framework completely different with respect to the 
traditional KS one. 

In this work we aim at extending our knowledge on 
the functional W a [p) by studying its behavior for a < 0, 
thus considering attractive electrons. In the very limit 
a — > — oo, we propose a rather simple and physically ap- 
pealing solution, which can be evaluated exactly in the 
case of two-electron systems. The main interest in the ex- 
ploration of this limit is to find new pieces of information 
on the unknown functional W Q [p]. Similarly to previous 
work on the opposite a — > oo limit, we expect that our 
results will open new possibilities for improving state-of- 
the art DFT approximations. Notice that, so far, the 
a < case has been only addressed in the special case of 
two electrons confined on the surface of a sphere, where 
accurate numerical calculations for a e (— oo, oo) have 
been performed [35]. A functional tuned to reproduce 
these calculations also for a < has been recently pro- 
posed [36]. The present work fully extends these early 
results by addressing the a — > — oo limit in a general 
way. 

The paper is organized as follows. In the next Sec. II 
we present our solution for ^ a [p] , t>" xt ( [p] , r ) and W a [p] in 
the a — > — oo limit, discussing the implications for DFT, 
as well as the limitations of our approach. As a first appli- 
cation, we show in Sec. Ill that for two-electron systems, 
the information on the a — > — oo limit can be used to 
get a rather accurate estimate of E2 h2 [p] of Eq. (10) by 
simply interpolating between the a — > ±oo limits, and 
thus without using the unoccupied orbitals. The same 
procedure is also applied to the Be isoclcctronic series, 
and the role of near-degeneracy effects in this framework 
is discussed. The last Sec. IV is devoted to conclusions 
and perspectives. 



II. THE LIMIT a -» -oo 

For a < 0, the Hamiltonian (7) describes attractive 
electrons with a given ground-state density p — p(r). 
To understand how attractive fermions can be forced to 
form a given smooth density distribution p by means of a 
local external potential v" xt ([p],r), we consider here the 
extreme case a — ► — oo. 
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FIG. 1: As the coupling strength a — > — oo, the N electrons 
are expected to form a compact (point-like) attractive elec- 
tron cluster (AEC), whose center of mass position has the 
probability distribution ip(r). 



A. The wave function f Q [p] for a — > — oo 

As their attraction becomes very strong (a — * — oo), 
we expect that the electrons in the state ^ a \fi] form a 
compact "attractive-electron cluster" (AEC). With this, 
we mean that simultaneous measurement of their posi- 
tions in this state will always yield N points ri,...,rjv 
in space that are very close to each other, much closer 
than any distance over which the density p(r) changes 
appreciably. For N — 2, e.g., the AEC is a positronium- 
type object with an average distance of (|ri — r 2 |) = j^- 

between the two electrons (where as = h 2 /m c e 2 is the 
Bohr radius). Notice that here we are only interested in 
the mathematical limit a — > — oo of the Hamiltonian (7), 
so that we disregard any relativistic effect. 

Exploiting this concept of a compact AEC, we expect 
that, as a — * — oo, the external potential w" xt ([p],r) ap- 
proaches a smooth function of r which gives the quasi 
point-like AEC the probability distribution -^p(r), as 
pictorially sketched in Fig. 1. In order to formalize this 
idea, we write 



N 

X] u Sct([p]> r i)*a([p];ri,...,rAr;o- 1 ,...,c7^) 
»=i 

~ ^«e X t([p]j R ) *a([p];ri, -,r N ;ai, ...,a N ) 

(a«-l), (15) 



where R = -k Yli=i v i IS the center of mass of the N 
electrons. The accuracy of this approximation grows in- 
definitely as a — > — oo, when the radius of the AEC tends 
to zero. Thus, we introduce relative Jacobi coordinates 



s„ (n = 1, N — 1) and the center-of-mass s^r = R 
»• - r n+ i (1 < n < N), 



s N 



l — l 

1 N 

t=i 



R. 



For 1 < n < N, the inverse transformation reads 

N-l 



R 



n - 1 



n 



E 



(16) 
(17) 

(18) 



while for n = 1 and n — N we have 

N-l 



e+i 



N - 1 

r N = R ^-sjv-i. (19) 



In terms of the Jacobi coordinates 8%, ...,Sjv_i,sjv = 



R, the operator f = — 52 i=1 Vf assumes the form 
h 2 d 2 



T 



T, 



rcl 



M = Nm e , (20) 



2M <9R 2 

where T re i acts on the relative coordinates only, 



2 ^ m n dsl ' 



n + 1 



-m„. 



(21) 



Also, the purely multiplicative operator V^o of Eq. (3) 
depends on the relative coordinates only, 



V cc = V cc (si, Sjy-l), 



(22) 



since — Tj is independent of R, see Eqs. (18) and (19). 

The resulting structure of the Hamiltonian (7) , within 
the approximation (15), implies a product ansatz for the 
wave function \l/ a [/>] in terms of the new coordinates, 



*a[p] 



i a ([/?]; R) tp a (si, -,Sjv-i; ai, a N ) (23) 



where <fi a and are, respectively, the lowest-eigenvalue 
solutions of the following Schrodinger Equations, 



• Nv. 



h 2 d 2 

2M&R 2 

|f rcl + aVee{si, ...,S N -i)\lp a 



SdMR)}< 



(24) 

• (25) 



Since 0Q,([p];R) is symmetric with respect to permuta- 
tions of the electronic coordinates r„ , the second factor 
"0q of the wave function (23) must be anti-symmetric, 

m-) = 7=r E (-irArva-)- (26) 

Here, Sat is the group of the AH permutations it of A*" 
elements, and (— l) ir is the sign of it. The operator is 
defined by 



iV0a:({s(ri, ...,riv)}; <7i, crjv) = 

V'a({s(r 7r ( 1 ), ...^^(jv))}; (7^(1), ov(jv)) 



(27) 
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with the short-hand notation {s} := si,...,sjv-i and 
the explicit functions s„ = s„(ri, rjv) of Eq. (16). 
To recover on the right-hand side a pure function of 
Si,...,sjv_i ( an d of the a n ), the r n r n \ must be re- 
expressed in terms of the s„ by Eqs. (18) and (19), where 
cancels out. 

The factors <f) a and ip a of the wave function (23) are 
addressed, respectively, in the following two subsections. 



B. The external potential v" xt ([p],r) for a — > — oo 

In terms of the first factor (f> a in Eq. (23), the proba- 
bility density of the electronic center-of-mass position R 
in the wave function ^f a [p] for a — > — oo is given by 

p£n([p],R) - J2 /<*V" jd^s^^M 2 

£7l,...,<7JV 

= l^([p],R)| 2 - (28) 

For a — > — oo, when the AEC becomes very compact 
(point-like), the function Np^ m ( [p] , R) must approach 
the electron density p(r). Then, using (f> a (r) = p" m (r), 
Eq. (24) can be resolved for v" xt ([p],r), 



lim 



N 



fi 2 V 2 7^) 



2m N 2 



(29) 



According to this result, the external potential w" xt ([/9], r) 
in the Hamiltonian (7) approaches a smooth and finite 
function v~^([p],r), as a — > — oo. The value -E<^° can 
be fixed by the condition w" xt ([p],r) — > for |r| — > oo. 
Obviously, Nv~^° ( [p] , r) is simply the external potential 
for a single particle with mass M — Nm c (which is the 
AEC of the previous Subsec. II A) and ground-state den- 
sity jfP(r). 

So far, we see that an external potential «ext([p]> r ) can > 
for a — > — oo, generate a given smooth density distribu- 
tion p(r) for N attractive electrons. The smooth distribu- 
tion is achieved by the uncertainty in the center-of-mass 
position R of the AEC. This picture becomes question- 
able in cases such as the stretched H2 molecule, where 
the density is separated spatially into several pieces. 



C. The integrand W a [p] for a — > — 00 

As we have seen in Subsec II A, as a — > —00 *f? a [p] 
should be better and better approximated by the prod- 
uct ansatz of Eq. (23). Consequently, the expectation 
(*a[p]|^ee|*a[p]) can be computed in this limit by using 
the second factor ip a of Eq. (23), as V ee only depends on 
the relative coordinates. Thus, for a — > —00 we have 

W a [p} + U[p] = (^ a [p}\V ee \^ a [p]) - 

J2 jd i S l --- y > d 3 S A r_ 1 |^ Q (...)| 2 Kc(si,...,S A r_ 1 ),(30) 



where we have integrated out the center of mass coordi- 
nate and used the normalization J d 3 R\<p a ([p], R)| 2 = 1. 
We now remark that the Schrodingcr equation for ip a 
of Eq. (25) is not affected by the details of the den- 
sity function p(r), but only by the electron number 
N = J d 3 rp(r), so that the same is true for its lowest- 
eigenvalue solution \p a and thus for the right-hand side 
of Eq. (30). This means that, for all different densi- 
ties p with the same electron number N, the function 
W a [p] + U[p] on the left-hand side of Eq. (30) has the 
same asymptote when a — > —00. This point is of course 
very appealing, but has obvious limitations that will be 
discussed in the next Subsec. II D. 

To address the solution ip a of Eq. (25), we consider the 
universal (i.e., density-independent) hamiltonian 



H" = T + aV ei 



(a < 0). 



(31) 



H£ describes a system of N attractive electrons in the 
absence of any localizing external potential v ext (r). Such 
a system has translational symmetry and a uniform 
ground-state density p. To make p finite, we can in- 
troduce a normalization volume Q — N/p and impose 
periodic boundary conditions on the wave function. In 
the ground state of , the N electrons are forming a 
"free" AEC whose center-of-mass position R has a uni- 
form probability distribution within £2, with density -^p. 

Notice that for the hamiltonian (31), Eq. (15) is not 
an approximation, since v" xt ([p],r) = 0. Consequently, 
the product on the right-hand side of Eq. (23) , where the 
factors 4> a and ip a satisfy Eqs. (24) [with v" xt = 0] and 
(25), respectively, represents the exact ground state 
of , for all finite values of a < 0. Therefore, Eq. (30) 
can be written as 

(* a [p]\V ee \* a [p])^{*a\Vee\^) (<*--«>). (32) 

Due to the universal form of the Hamiltonian (31), its 
ground-state wave functions for different interaction 
strengths a < are related by a simple scaling law, 

*»({r,a}) = \ar^_ 1 ({\a\r,a}), (33) 

in the short-hand notation {r, a} = ri, rjy; <7i, oat. 
The resulting scaling law for the ground-state energy E£ 
of H£ reads 

^ = « 2 ^=-l- (34) 
Notice that E£ is the eigenvalue in Eq. (25), 

{f re i + aVee(8l,-,8iV-l)}Va = ^aV'a ■ (35) 



The virial theorem, (tf£|T|tf£) = -i(*^|at> e e|*a ). 
yields E% = i<*£|avg*£> or 



(n\v~\n) = -K = **E% 

a 



(36) 
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Consequently, according to Eq. (32), the integrand of 
Eq. (1) asymptotically approaches for a — > — oo a lin- 
ear function with slope 2E^ 1 < 0, 



W a [p] 



2a E 



N 



U[p] (a — ► — oo). 



(37) 



As said, since E T ^ 1 is the ground-state energy of a free 
AEC at a = — 1, it cannot depend on details of the den- 
sity p, but only on the electron number N = J d 3 r p(r). 
We discuss below the limitations of this result. 



D. Limitations of the AEC solution 

First of all, as already mentioned, the AEC picture 
breaks down when the density is composed of several 
disjointed subsystems (e.g., stretched molecules). In the 
case of non-interacting fragments, in order to be size con- 
sistent, the limit of Eq. (30) should apply separately to 
each subsystem, which is physically reasonable, but dif- 
ficult to prove rigorously. The transition from a jointed 
to a disjointed system (e.g., during molecular stretching) 
remains also unclear. 

An even more important point is the following. In or- 
der to be useful to construct approximations for 
the a — ► — oo asymptote of Eq. (30) should be approached 
in a continuous and smooth way starting from a = 1 (i.e., 
without phase transitions). This is generally not true, so 
that often the state ip a which is smoothly connected with 
the physical (a = 1) system is not the ground state of 
the hamiltonian of Eq. (25). As a simple example, con- 
sider the N atom. At a = 1, we have JV-j = 5 spin-up 
electrons and TVj = 2 spin-down electrons, but we may 
expect that the ground state of the attractive electron 



cluster corresponds to N-t 



4, N 



3, as pairing is 



usually energetically advantageous for an attractive in- 
teraction, so that we will have a spin-flip transition at 
a < 0. This simple case might be solved in a spin-DFT 
framework (in which we look for the a < systems with 
the same spin densities of the physical system) , but we 
may encounter much more difficult cases as the number 
of electrons increases, with possible "exotic" phases. An 
interesting point is that, however, the a — > — oo potential 
v" xt ([p], r) of Eq. (29) is the same for all the excited states 
of the cluster, including the one continuously connected 
with the physical system, if such a state exists. Nonethe- 
less, with all these limitations in mind, we may still ex- 
pect that the results presented in the previous section be 
qualitatively correct for many interesting systems. 

In simple cases such as closed-shell two- and four- 
electron systems with a reasonable compact density, we 
might expect a continuous connection between the limit 
of Eq. (37) and the a = 1 case. In the next Sec. Ill, we 
construct a simple continuous interpolation between the 
two limits a — > ±oo and we use it to estimate E^ L2 [p] for 
spherical two- and four-electron densities. In the latter 
case, we will discuss the role of near-degeneracy effects. 



III. £c [p] FROM THE LIMITS a -> ±oo 

We now construct an analytic model VF AR [p] ("AR" 
stands for "attraction-repulsion" here) that shares the 
properties (13), (14), and (37) with the unknown exact 
integrand W a [p) of Eq. (1), which we assume here to be 
a continuous and smooth function of a also for a < 0, 
even if we keep in mind the limitations of this assumption 
discussed in the previous Subsec. II D. In particular, our 
model will share the correct density-scaling behavior with 
the unknown exact integrand W a [p], see Eq. (53) below. 
Since the inverse function is easier to model, we write 



W^[p] = f-\a). 



(38) 



According to Eqs. (13), (14), and (37), respectively, this 
inverse function /(to) must satisfy the conditions 



fH 

f(E x )= 

/M 



WL 



U 



2E 



N 



(w^Woo), (39) 
(w = E x ), (40) 
(w — > oo). (41) 



Here, we have dropped the functional symbol [p] for 
brevity, as we shall do in most of the following equations. 

Eqs. (39) and (41) are satisfied by the preliminary 
choice 



fiH = 



w + U 
~2E W ~ 1 



+ 



W 



w-W a 



(42) 



which, however, ignores condition (40). In particular, 
fi(w) approaches the asymptotic form of Eq. (41) too 
slowly, fi(w) = (to + U)/2E N 1 + O(to~ 2 ). Instead, we 
expect /(to) = (to + U)/2E N 1 + 0(e~ w ) or, equivalently, 



W a [p] 



2a E" 



U[p] + 0(e a ) (a-»-oo). (43) 



This conjecture is based on Eq. (33), indicating that the 
size of the AEC shrinks by the factor |a| _1 as a — > 
— oo, while the external potential u" xt ([ j o],r) approaches 
a smooth finite limit, sec Eq. (29). Consequently, the 
point where the radius of the AEC becomes smaller than 
any distance |r| over which v" xt ([/?], r) changes apprecia- 
bly, must be reached at some moderately negative value 
ao[p] of the parameter a, 



W a [p] w 2aE^-U[p] (a<a [p)). 



(44) 



This conjecture is also supported by the numerical data 
sets shown in Fig. 3, where ao takes values w —3.4 ('Sp3') 
and w -1.0 ('Sp2'). Accordingly, we modify Eq. (42), 



/2H 



w + U 
~2E^ 



W 



w-W a 



e-( w ~ w ~y r . (45) 



Note that the exponential factor in Eq. (45) equals 1 at 
w = Woo[p], thus preserving the correct limit (39). Due 
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FIG. 2: The functional W a [p] as a function of a (solid curve), 
as modeled in Sec. Ill, using the five parameters for the Be 
atom in Table I. 



to the general scaling property of the exact W a [/t)], see 
Eq. (53) below, the parameter r — r[p] must scale the 
same way as Woo[p]j F[p\] — \r[p]. Therefore, we set 



r[ P ] = | Woo [p]|, 



(46) 



identifying the unknown number ao[p] from Eq. (44) ap- 
proximately with the zero of W a [p] , 



W ao[p] [p] « 



/(0) « a [p], (47) 



see Figs. 2 and 3. 

To satisfy condition (40), too, we introduce an extra 
factor hs(w) with /is(Woo) = 1, 



, . w + U 



W' 

Y V OQ 

w-W a 



In terms of the number 
E x + U /R 



h B (w) e -(™-w-)/r i 
(48) 



B = 



-2E N 

we choose for w 6 [Woo , +oo] 



L ( E * W W °° ) e ^- w ^ r > 0, (49) 

1 V "oo / 



h B (w) 



l + (5-l)fE^ 

w-W x 



B > 1, 
B < 1. 



(50) 



Then, /i B (Woo) = li M^*) = S, and h B {w) > is 
monotonically increasing (B > 1) or decreasing (B < 1). 

The function W^ R [p] = / _1 (a) has five parameters: 
E N _ l , U[p], E x [p], Woolfi], and W^J/o], and it is shown in 
Fig. 2. As an illustration, Fig. 3 shows the true integrand 



FIG. 3: The functional W^p] as a function of a for the sys- 
tems 'Sp2' and 'Sp3' (see Appendix A). The accurate nu- 
merical values (dots) are compared with the model of Sec. Ill 
(solid curves). 



VFcJp] (dots), evaluated numerically for the systems 'Sp2' 
and 'Sp3' (see Appendix A). For these systems, the 
evaluation of W Q [/o] is straightforward, since their (in- 
dependent external potential (r) in the Hamiltonian 
(7) is trivial and thus known from the beginning. 

Since 2£ c GL2 [p] = £W a [p]\ a=0 « !//'(£*)> our model 



W^ R [/o] yields a simple prediction E^ L2 [p] = 2 ^ E ) 
the second-order correlation energy E2 L2 [p], 



for 



E^[p] = E 



N 



1 + (E x + U) 



B' 



B' 



-Woo 

3-B 



1 

" Woe 

(B>1), 
(B<1). 



(51) 



Unlike the exact functional E^ L2 [p] whose evaluation re- 
quires additional knowledge of all the unoccupied KS or- 
bitals, the approximation (51) only depends on the N 
occupied orbitals (via E K [p]), on the universal number 
E_±, and explicitly on the density p itself. Nevertheless, 
the new functional has the correct scaling behavior of the 
exact E2 L2 [p] (only valid when the KS system does not 
become degenerate or quasi-degenerate), 



E^ 2 [p x ] = E™[p] 



(52) 



where p\(r) = A 3 p(Ar) (with A > 0) is a scaled density. 
More generally, the model integrand \p] has the cor- 
rect scaling property [37] of its exact counterpart W a [p] , 



W^[ PX ) = \W$\p] 



(53) 



This law has a simple graphical interpretation: Plotting 
[p\] versus a amounts to zooming the corresponding 
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plot of VK^ R [p] by the factor A. This graphical property 
is satisfied by the function f(w), Eq. (48) (and then also 
for its inverse, W^ R [p]), since we have the general scaling 
laws U[p x ] = \U[p], E x [p x ] = \EM, W 00 [p x ]=XW 00 [p], 
and V^M = A 3 / 2 F<[p]. 



A. Application to N = 2 systems 

For a closed-shell two-electron system with a reason- 
able compact density the cluster energy E^ 2 needed 
in Eq. (37) can be calculated exactly: the AEC defines 
a positronium-like problem whose ground state energy is 
E 2 _ 1 = j Ha. Moreover, in this special case we expect 
that the a — > — oo limit is reached smoothly from a = 1. 

In Table I we report the prediction of E^ h2 [p] from 
Eq. (51) for several two-electron densities, and we com- 
pare it with accurate values from the literature. The 
systems denoted 'Sp2' and 'Sp3' correspond to electrons 
confined on the surface of a sphere (in 2D and 3D, respec- 
tively) and are defined in Appendix A. For these systems 
exact values are available [35, 38, 39]. 'Exp' refers to 
a hypothetic two-electron atom with ground state den- 
sity p(r) = e~ r / aB /4na%, whose corresponding accurate 
E^ h2 [p] is taken from Ref. 40. 'Hoo' refers to the Hooke 
atom, consisting of two electrons in the external potential 
^ext( r ) = §?* 2 > w hh k = \. Its exact density is analyti- 
cally known [41], and the corresponding E^ L2 [p] is from 
Ref. 42. The accurate densities for He and Ne 8+ and the 
corresponding values of E^p L2 [p] are taken from the work 
of Colonna and Savin [43]. For all the densities consid- 
ered here, the junctionals Woo[p] and W^p] have been 
evaluated using the procedure described, respectively, in 
Refs. 26 and 27. 

We see that for all the two-electron systems consid- 
ered here, the functional E2 L2 [p] of Eq. (51) is in very 
good agreement with the corresponding accurate values 
from the literature, with a maximum error of ~ 3 mHa. 
Notice in particular how, for the 'Sp2' system, the un- 
usually high value E^ L2 [p] = —227 .4 mHa of the slope 
i^W a [p]| a= o in Fig. 3 is predicted from the unusually 
high energy E 2 _ 1 = — 1 Ha (i.e., the slope of the corre- 
sponding inclined dashed line at a < in Fig. 3) of the 
AEC in two dimensions. In contrast, the remaining four 
parameters U, E x , Woo, and in Table I do not differ 
dramatically between 'Sp2' and 'Sp3'. 



B. Application to the Be series: 
The effect of near-degeneracy 

For the Be isoelectronic series, we have closed-shell 
N = 4 systems, so that we may still expect that the 
a — > — oo limit of Eq. (37) is approached smoothly from 
a = 1. However, as the nuclear charge Ze increases, the 
2s and 2p KS orbital energies become closer and closer, 
resulting in a large increase of E^ L2 [p\. Indeed, in this 
case the scaling property of Eq. (52), satisfied by the 



present model of Eq. (51), does not hold anymore, so 
that we can expect a failure of Eq. (51) for large Ze. 

As an illustration, we have applied the model of 
Eq. (51) to the Be and the Ne 6+ densities. The cluster 
energy for N = 4 has been estimated from a configuration 
interaction (CI) calculation [44]. The accurate values of 
U[p], E x [p] and E2 h2 [p] are taken, in both cases, from 
the work of Colonna and Savin [43] . Using the same den- 
sities of Colonna and Savin [43], we have calculated here 
the values of [p] and [p] , following the procedures 
described, respectively, in Refs. 26 and 27. 

We see that the result for Be is quite accurate, yield- 
ing an estimate of E^ h2 [p] with an error of only 2 mH, 
confirming the hypothesis that, also in this case, the 
a — > — oo limit is reached smoothly from a = 1. However, 
as the system starts to display strong near-degeneracy 
as in the case of Ne 6+ , we see that the present model 
fails, making an error of 60% in the estimate of ££* L2 [p]. 
When Ze continues to increase beyond the value Z = 10 
considered here, the ground-state density becomes no 
longer pure-state w-representablc, but only ensemble v- 
representable [45-50]. At this point, the definition of 
both E x [p] and E^ 2 [p] should change. 



IV. CONCLUSIONS AND PERSPECTIVES 

We have presented a comprehensive analysis of the 
DFT adiabatic connection at negative coupling strengths 
a. In the extreme limit a — ► — oo we have found a 
simple and physically appealing solution, which, albeit 
suffering from limitations due to the complexity of the 
many-electron problem with attractive interaction, can 
be calculated exactly in the case of reasonably com- 
pact two-electron systems, and provides insight also when 
N > 2. As a first example of application of our results, 
we have shown that, for N = 2, the exact information at 
a — ► — oo can be used, in combination with the opposite 
a — > +oo limit, to estimate the second-order correla- 
tion energy E^ 2 [p] without using virtual orbitals. The 
same procedure works very well also for the Be atom, 
but breaks down in the case of Ne 6+ , because of strong 
near-degeneracy effects. 

The analysis carried out here extends our knowledge 
on the exact properties of the exchange-correlation func- 
tional E^c[p], and can be used to test and improve 
approximations. Another possible application of the 
present results could arise in the framework of a recently 
proposed approach to the many-electron problem that 
contains a mixture of Hartree-Fock and Hartee-Fock- 
Bogoliubov methods [51]. This approach is based on 
the splitting of the Coulomb electron-electron interac- 
tion 1/r as — l/r+2/r. The attractive part — 1/r triggers 
the Hartee-Fock-Bogoliubov solution. An energy density 
functional to describe dynamical correlation is also added 
to the method [51]. For this latter point, the extension 
to negative coupling a of the adiabatic connection could 
prove useful in defining and constructing the appropriate 
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System 


N 


D 


E N , 


U[ P ] 


EM 




B 


£c GL2 [p] 


£ C GL2 [P] 


error 


Units 


- 


- 


1 Ha=27.21 eV 


- 


lQ-'- 5 Ha 


10 _a Ha 


'Sp2' 


2 


2 


-1.00 


2.000 


-1.000 


-1.500 


0.250 


2.791 


-228.1 


-227.4 


-0.7 (0.3%) 


'Sp3' 


2 


3 


-0.25 


1.698 


-0.849 


-1.198 


0.375 


1.968 


-46.5 


-47.6 


+1.1 (2.3%) 


'Exp' 


2 


3 


-0.25 


1.250 


-0.625 


-0.910 


0.345 


1.167 


-43.4 


-46.7 


+3.3 (7.1%) 


'Hoo' 


2 


3 


-0.25 


1.030 


-0.515 


-0.743 


0.208 


1.682 


-47.2 


-50.5 


+3.3 (6.5%) 


He 


2 


3 


-0.25 


2.049 


-1.025 


-1.500 


0.621 


1.649 


-48.6 


-47.5 


-1.1 (2.3%) 


Ne 8+ 


2 


3 


-0.25 


12.055 


-6.028 


-8.794 


8.792 


1.631 


-48.0 


-46.7 


-1.3 (2.8%) 


Be 


4 


3 


-1.255 


7.217 


-2.673 


-4.021 


2.59 


0.6857 


-126.4 


-128.4 


+2.0 (1.6%) 


Ne 6 + 


I 


3 


-1.255 21.742 


-7.600 


-11.563 


12 


0.8655 


-127.5 


-320.4 


+193 (60%) 



TABLE I: The AEC energies and the coefficients U[p], E x [p], W x [p], and W^[p] of the expansions (37), (13), and (14), 
the dimensionless parameter B, Eq. (49), and the resulting predictions E^ h2 [p] of the second-order correlation energy E^ h2 [p] 
for various D-dimensional iV-electron systems. The systems 'Sp2' and 'Sp3', here with radius R — as, are defined in Appendix 
A; 'Exp' refers to a hypothetic two-electron atom with ground state density p(r) = e~ r l aB /4:Tva 3 B ; 'Hoo' refers to the Hooke 
atom, consisting of two electrons in the external potential w cx t(r) = |r 2 , with k = \. 



energy functional needed in the method. 



with R = as in Table I. 
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APPENDIX A: TWO ELECTRONS ON THE 
SURFACE OF A SPHERE 



A similar, but more realistic two-electron system is 
'Sp3' which has its electrons confined to the 3D surface 
1S3 of a sphere in 4D space, 



S 3 = {(x,y,z,u) e R i \x 2 + y 2 + z 2 + u 2 = R 2 Y (A4) 



'Sp2' is a system of two electrons that are confined 
to the 2D surface 1S2 of a sphere with given radius R. 
Recently [35, 38, 39], this system has been studied ex- 
tensively. It is of particular interest for DFT, since its 
2D ground-state density for all values of the interaction 
strength a £ R is distributed uniformly over the spherical 
surface. Consequently, the external potential v" xt ([/j];r) 
in the Hamiltonian (7) is trivial and known from the be- 
ginning. The 2D confinement of these electrons, however, 
gives rise to an unusually large value [35, 38] of 

£ c GL2 (i?) = (4 In 2 - 3) Ha = -227.4 mHa. (Al) 

This effect can be traced (see Fig. 3) to the strongly neg- 
ative energy of a 2D AEC [35] , 



E 



N=2 



1 Ha (2D), 
(3D). 



-| Ha 



(A2) 



While the value (Al) is independent of the spherical ra- 
dius R, the remaining functional of Table I are 



E X (R) 

wUR) 



~R = 

2 R' 
1 /a B \ 3 / 2 e 2 
4 V ~R~) ~ob 



U(R), 



Again, the ground-state density is distributed uniformly, 
but now in a 3D finite (curved) space S3. Consequently, 
we obtain a value very close the one for the He atom (see 
Table I) [39, 52] 



E? L2 (R) = -47.6 mHa. 



The remaining functionals of Table I are [52] , 



E X (R) 
W^R) 

WL(R) 



j^e? _1 
3tt R " 2 

1 _ 16\ e 2 

2 3n) H 

e 

<>B 



U(R), 



3 /a B \3/2 e 2 
8 V^R - 



(A3) with R = a B in Table I. 



(A5) 



(A6) 
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